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2012-2013  Annual  Progress  Report 

This  report  presents  the  specific  aims  and  accomplishments  of  our  prostate  cancer  research  project  during 
the  year  of  funding  sponsored  by  the  US  Department  of  the  Army.  It  covers  our  activities  from  May  2012  to 
June  2013. 

1.  Summary  of  our  Year  1  (2011-2012)  research  and  development: 

We  have  designed,  implemented,  and  tested  a  frequency-domain  photoacoustic  imaging  system  using 
tissue-like  phantoms.  The  imaging  principle  and  data  processing  method  have  been  presented  and 
verified.  The  measured  photoacoustic  (PA)  signal  strength  depends  on  the  absorption  coefficient,  depth 
of  the  target,  and  the  modulated  light  intensity  of  the  laser  diode.  We  have  learned  that  appropriately 
processing  the  measured  PA  data  and  background  interference  is  an  important  key  to  correctly  display 
the  PA  images.  The  lateral  spatial  resolution  of  PA  images  depends  on  ultrasound  frequencies  and 
focal  size  of  the  ultrasound  transducers.  Although  the  current  system  has  a  low  axial  resolution  (due  to 
the  adopted  ultrasound  transducer),  the  axial  resolution  can  be  improved  by  using  a  high  numerical 
aperture  ultrasound  transducer  [1]  or  by  adopting  a  frequency-swept  technique  [2].  Our  results  and 
conclusion  are  accepted  for  a  peer-reviewed  publication  [1],  which  is  attached  as  an  attachment. 

In  the  meantime,  however,  we  also  learned  that  the  frequency-domain  photoacoustic  imaging 
method  developed  by  our  approach  has  a  severe  limitation  in  sensing  tissue  signals  deeper  than  5  mm. 
Then,  we  tried  to  alter  or  modify  our  original  idea  by  developing  an  integrated  approach  between 
ultrasound  and  optical  tomography,  namely,  transrectal  ultrasound-guided  diffuse  optical  tomography 
(TRUS-DOT)  for  improved  prostate  cancer  imaging.  In  the  last  10-12  months,  we  have  made  great 
progress  and  obtained  promising  results  using  tissue-phantom  experiments.  We  believe  that  our 
current  development  can  lead  to  another  innovative  imaging  method,  which  can  correctly  detect  or 
identify  aggressive  prostate  cancer  tissue  from  non-aggressive  tissues.  In  this  report,  we  provide  details 
on  this  modified/updated  approach  with  computer  simulations  and  experimental  results. 

2.  Body  of  the  Report 

2.1  Hierarchical  clustering  to  improve  TRUS-DOT  for  prostate  cancer  imaging 

Diffuse  optical  tomography  (DOT)  is  a  non-invasive  imaging  modality,  which  utilizes  near  infrared  (NIR) 
light  and  provides  blood-based  absorption  maps  within  the  detected  tissue  volume  interrogated  by  the 
optical  optodes.  When  using  multiple  wavelengths,  DOT  is  also  capable  of  measuring  chromosphere 
concentrations,  such  as  oxy-hemoglobin,  deoxy-hemoglobin,  and  water.  Usage  of  DOT  for  breast 
cancer  detection  has  been  extensively  studied.  DOT  instrumentation  can  be  divided  into  three 
categories  based  on  the  principle  of  operation:  (1)  time  resolved  systems,  (2)  frequency  domain 
systems  and  (3)  continuous  wave  (CW)  systems.  Measurements  in  DOT  can  be  made  in  the  form  of 
transmission,  reflectance,  or  both.  Time  resolved  systems  rely  on  photon  counting  or  gated  imaging, 
which  measures  time  of  flight  of  photons  through  the  tissue.  However,  these  systems  are  very  costly 
and  bulky  in  comparison  with  CW  systems.  The  frequency  domain  systems  modulate  the  light  source 
(e.g.,  laser  diode)  typically  in  the  radio  frequency  range  (e.g.,  100  MHz)  and  measure  the  amplitude 
and  phase  shift  of  the  detected  signals.  CW  systems  are  the  simplest,  cost  effective,  and  fastest  in 
data  collection;  they  can  be  made  to  image  tissues  at  a  video-rate.  However,  since  CW  systems 
measure  only  intensities  of  the  reflected  or  transmitted  light,  they  cannot  separate  the  absorption  and 
scattering  effects  of  the  tissue. 

Despite  of  the  instrumentation  utilized,  DOT  has  poor  spatial  resolution.  Measurements  in 
reflectance  geometry  are  even  worse  in  spatial  resolution  than  in  the  transmission  geometry.  The  cause 
that  limits  the  spatial  resolution  of  DOT  is  light  scattering  and  diffusion  in  tissues,  making  DOT  have  to 
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solve  an  underdetermined  and  ill-posed  inverse  problem.  One  way  to  improve  the  spatial  resolution  is 
to  couple  DOT  techniques  with  other  imaging  techniques,  such  as  MRI,  ultrasound  (US),  and  others. 
This  process  is  usually  done  by  combining  a  priori  information  in  the  inversion  procedure.  In  particular, 
a  combined  usage  of  TRUS  with  DOT  could  improve  the  accuracy  in  identifying  prostate  cancer  1-2 
centimeters  below  the  prostate  surface,  as  studied  previously  [3],  However,  this  reported  approach 
treated  each  anatomical  region,  such  as  the  rectum  wall,  prostate  gland,  and  prostate  cancer  lesion,  to 
be  homogenous  and  thus  to  be  reconstructed  with  uniform  optical  properties  within  each  region. 
However,  the  reconstructed  images  by  the  reported  method  [3]  could  be  erroneous  because  of  the  fact 
that  TRUS  does  not  provide  a  clear  or  clean  signature  for  prostate  cancer.  On  the  other  hand,  it  has 
been  demonstrated  that  the  spatial  resolution  of  reconstructed  DOT  images  depends  upon  the  number 
and  locations  of  optical  sources  and  detectors  utilized.  Given  the  limited  space  available  in  TRUS,  the 
idea  to  incorporate  more  sources  and  detectors  for  a  combined  TRUS-DOT  approach  seems  not  to  be 
very  practical. 

To  solve  the  problems  given  above,  we  made  new  development  in  two  aspects:  (1)  development  of 
hierarchical  clustering  method  (HCM)  for  improved  contrast  recovery  and  spatial  resolution  in  DOT  for 
prostate  cancer  imaging;  (2)  development  of  Trans-rectal  compatible  optical  probe  and  instrumentation 
for  prostate  cancer  imaging.  First,  in  Section  2.2,  we  will  explain  the  HCM  procedure  and  demonstrate 
the  capability  of  HCM  using  computer  simulations.  Then,  in  Section  2.3,  we  will  explain  the 
instrumentation  utilized  for  trans-rectal  compatible  DOT  and  demonstrate  the  capability  of  HCM  using 
laboratory  phantoms. 

2.2  Introduction,  implementation,  and  Investigation  of  hierarchical  clustering 

Principle  and  algorithm  of  HCM 

In  standard  practice,  the  equation  used  to  iteratively  update  DOT  image  reconstruction  is  given  by  [4,5] 

(JTJ  +  2AI)(cpa)=JT(y-F(Ma))  (1) 

where  J  is  the  Jacobin  matrix  (or  called  sensitivity  matrix),  /  is  the  identity  matrix,  y  is  the  measured 
data  in  a  matrix  form,  F  is  the  operator  that  generates  the  forward  model  for  light  propagation  in  tissues, 
pa  is  the  vector  of  optical  properties,  and  A  is  the  regularization  parameter.  Note  that  changes  only  in  pa 
are  considered  in  eq.  (1),  because  our  DOT  measurement  utilizes  CW  NIR  light  with  an  assumption 
that  variation  in  light  scattering  across  the  medium  is  minimal. 


In  our  HCM  method,  we  divide  the  region  of  interest  (ROI)  (i.e. ,  the  prostate  region)  into  several 
geometric  clusters.  We  assume  that  each  of  the  geometric  clusters/units  is  homogeneous  and  has  the 
same  optical  property.  In  this  way,  the  computational  domain  could  be  partially  heterogeneous  since 
the  domain  may  contain  several  geometric  units.  Specifically,  the  nodes  in  the  mesh  are  tagged  and 
separated  into  subsets  (GSi,  GS2...GS/')  with  respect  to  each  unit.  The  Jacobian  matrix  in  eq.  (1)  is 
then  modified  to  J*  and  is  given  by  where  S  matrix  had  a  size  of  NN*NS,  where  NN  is  the 

number  of  nodes  within  the  specific  unit/cluster  and  NS  is  the  number  of  geometric  clusters.  The 
elements  of  S  matrix  are  given  by  [5]: 


HiJ) 


1  if  is  GS, 
0  else 


(2) 


At  the  end  of  each  iteration,  the  solution  vector  is  mapped  back  to  each  node  using  the  following 
equation 


8pa=S{8pa) 


(3) 
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/va*  is  a  vector  with  optical  properties  of  geometric  clusters.  When  reconstructing  two  ROIs  or  clusters, 
such  as  background  and  anomaly,  the  background  mesh  was  geometrically  segmented  in  a 
heterogeneous  fashion.  For  multiple  clusters,  we  hierarchically  implemented  the  proposed  method  by 
segmenting  the  region  which  is  more  prone  to  cancer,  utilizing  available  prior  spatial  information. 
Specifically,  we  applied  the  proposed  method  in  four  steps,  as  shown  in  Fig.  1. 


Reconstruct  optical  properties  for 
available  regions  (initial  guess  for 
step  2  and  step  3 


Divide  the  mesh  into  geometric 
segments 


Average  all  16  images  and  use 
FWHM  to  identify  suspicious 
gments 


reconstruct  -* 


Group  all  non  suspicious  segments 
as  region  and  suspicious  into 
_ several  regions _ 


Average  all  images  and  use 
FWHM  to  identify  suspicious 
segments 

Use  the  location  of  the 
suspicious  region  and  apply 
hard  prior  method 


reconstruct 


Stepl 
Step  2 
Step  3 
Step  4 


Figure  1 :  Flow  chart  indicating  steps  involved  in  HCM. 


In  step  1,  the  reconstruction  is  done  based  on  TRUS  images  and  the  assumption  of  a 
homogeneous  prostate.  With  such  prior  spatial  information,  the  reconstructed  iJa  values  in  both 
background  and  prostate  region  should  be  reasonably  accurate  with  respect  to  the  actual  values.  Then, 
the  reconstructed  /ja  values  in  available  ROIs  would  serve  as  the  initial  guess  in  step  2  and  step  3. 

Step  2  of  our  method  is  dedicated  to  finding  the  probable  locations  of  anomalies  (i.e.,  prostate 
tumors).  To  achieve  this,  we  geometrically  divide  the  prostate  region  into  several  clusters,  so  the 
prostate  tissue  becomes  a  heterogeneous  medium.  However,  without  prior  knowledge  of  suspicious 
locations,  dividing  the  prostate  into  several  clusters  may  result  in  mixing  of  suspicious  tissue  with 
normal  prostate  tissue  and  vice  versa.  In  order  to  prevent  this  situation,  we  segment  the  prostate 
volume  with  different  unit  volume  sizes  iteratively.  Specifically,  we  chose  a  unit  volume  as  an  initial 
segment  size  which  is  1x1x1  cm"  in  each  of  the  x,  y  and  z  dimensions.  Then,  we  vary  the  volume  size 
by  increasing  the  linear  length  of  the  volume  unit  in  each  dimension  iteratively.  For  example,  an 
increase  in  length  of  0.5  cm  in  only  x  direction  gave  rise  to  a  unit  volume  of  1 .5  x  1  x  1  cm3,  followed  by 
the  same  length  increase  in  only  y  or  only  z  direction.  In  this  way,  we  can  generate  8  different  unit 
volumes  in  three  x,  y,  z  directions,  by  increasing  the  linear  length  in  only  one  dimension  (x,  y,  z),  or  in 
two  dimensions  (  xy,  yz,  xz),  or  in  three  dimensions  (xyz)  and  the  initial  size. 

The  procedure  is  given  as  follows  (see  Fig.  1):  (1)  reconstruct  an  initial  \ja  image  with  a  starting 
base  unit  size  (i.e.,  1x1x1  cm3),  (2)  save  the  reconstructed  image,  and  go  back  and  change  the  unit 
volume  size  (e.g.,  1. 5x1x1  cm3  or  1.5x1 .5x1  cm3  or  1.5x1 .5x1 .5  cm3)  and  reconstruct  the  image  again 
(Step  2  in  Fig.1 ).  To  be  more  comprehensive,  we  next  increase  our  base  unit  volume  from  1x1x1  cm3  to 
1.5x1 .5x1 .5  cm3  and  then  to  2x2x2  cm3,  with  the  same  length  interval  of  0.5  cm  applied  to  increase  the 
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base  unit  (e.g.,  2x1 .5x1 .5  cm3,  2x2x1 .5  cm3,  or  2x2x2  cm3).  In  this  way,  we  are  able  to  generate 
another  set  of  8  reconstructed  images  with  varied  base  unit  volumes,  resulting  in  an  overall  16  (i=16) 
images  by  the  end  of  step  two.  (3)  We  average  all  of  the  reconstructed  \ja  images  to  obtain  the  final 
image.  (4)  Next,  we  search  for  suspicious  clusters  using  full  width  half  maximum  (FWHM)  of  the 
updated  pa  values  within  the  prostate  region.  If  no  suspicious  segment  is  identified,  we  conclude  that 
the  prostate  has  a  low  probability  of  having  cancer. 

In  step  3,  if  some  suspicious  clusters  in  step  2  are  seen,  we  then  group  all  the  non-suspicious 
clusters  as  one  new  single  segment  and  subdivide  the  suspicious  clusters  into  further  smaller  clusters. 
Now,  the  initial  unit  volume  size  used  within  the  suspicious  regions  in  step  3  is  set  to  0.5  x  0.5  x  0.5 
cm3.  The  procedure  explained  in  step  2  is  repeated  here  with  a  length  variation  of  0.25  cm  in  any  one  of 
three  dimensions.  Similar  to  step  2,  the  final  reconstructed  image  of  step  3  will  be  an  average  of  the  8 
images  (j=8)  obtained  by  varying  the  unit  volume  in  8  different  fashions.  FWHM  of  the  pa  values  is  still 
used  to  localize  suspicious  regions  for  further  inspection  with  an  improved  spatial  resolution. 

Finally  in  step  (4),  the  hard  prior  method  was  used  to  improve  the  optical  properties  of  the 
suspicious  regions  identified  in  step  3.  The  initial  guess  here  is  the  uniform  initial  guess  as  in  step  1. 

Computer  simulations 

To  validate  HCM,  we  performed  computer  simulations  by  considering  a  mesh,  which  was  anatomically 
similar  to  a  TRUS  prostate  image,  consisting  of  four  different  ROIs,  such  as  prostate,  peri-prostate 
tissue,  rectum  wall  and  prostate  tumors.  Two  spherical  shaped  inclusions  were  placed  with  a  horizontal 
separation  of  20  mm  and  a  vertical  depth  of  20  mm  from  the  surface  of  the  rectal  wall  (Fig.  2).  The 
rectum  wall  had  a  thickness  of  5  mm  and  a  curvature  radius  of  50  mm. 


Figure  2:  (a)  Probe  geometry  used  in  the  simulation;  each  optode  is  bifurcated  to  serve 
as  a  source  and  detector,  (b)  Two  anomalies  separated  by  20  mm  at  a  depth  of  20  mm. 
(c)  Mesh  has  been  sliced  isometric  view  to  show  the  simulation  geometry. 


In  step  1,  as  the  prior  location  of  prostate  cancer  is  often  not  visible  in  TRUS  images,  we  assumed 
homogeneous  optical  properties  for  all  ROIs.  A  simulated  TRUS-DOT  probe  was  used,  having  16 
bifurcated  optodes  that  could  serve  as  both  sources  and  detectors  (Fig.  2a).  The  CW  mode  was  utilized 
in  the  simulations;  1%  random  noise  was  added  to  the  data  to  mimic  the  instrument  noise.  Simulated 
NIR  data  was  computed  along  the  rectum  boundary  using  a  FEM-based  diffusion  forward  model  with 
NIRFAST  [6],  The  first-step  reconstruction  was  done  using  the  initial  guess  of  pa=0.01  mm'1  to  recover 
the  optical  property  of  the  prostate  (Fig.  3a).  Then,  the  reconstructed  data  from  the  first  step  served  as 
the  initial  guess  for  the  second  and  third  steps,  as  given  in  Section  2.  All  reconstructions  would  stop 
when  the  change  in  projection  errors  is  less  than  2%  of  the  previous  iteration.  Next,  the  prostate  tissue 
was  divided  into  several  clusters.  As  noted  earlier,  the  final  image  obtained  after  step  2  was  the 
average  over  16  independently  reconstructed  images,  as  demonstrated  in  Fig.  3b.  In  this  simulation, 
we  found  an  estimated  suspicious  region,  as  marked  by  the  dotted  rectangle  in  Fig.  3b,  using  FWHM. 
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In  step  3,  the  suspicious  region  identified  in  step  2  was  further  divided  into  smaller  clusters.  The  final 
image  of  step  3  shown  in  Fig.  3c  was  the  average  over  8  reconstructed  images.  FWHM  was  used  and 
resulted  in  two  suspicious  regions,  which  were  treated  as  two  individual  regions  and  entered  in  the  hard 
prior  method  for  finer  reconstruction  in  step  4. 


mt:  . 

'  '  ' _ ' 

r  ; 

: 

(d) 

Figure  3:  reconstructed  pa  values  in  mm'1  using  different  reconstruction  steps.  The  dotted 
circles  indicate  the  real  locations  of  the  anomaly.  Reconstructed  image  (a)  after  step  1  using 
HCM,  (b)  after  step  2  using  HCM;  dotted  rectangle  indicates  the  suspicious  location  selected 
using  FWHM,  (c)  after  step  3  using  HCM,  (d)  after  step  4  using  HCM.  (e)  Reconstructed  image 
without  any  inclusions,  but  reconstructed  using  HCM.  (f)  Reconstructed  image  for  the  same 
case  using  a  known  hard  prior  for  the  inclusions. 


2.3  Measurement  system  and  data  processing 
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Figure  4:  Instrumentation  and  probe  setup  utilized  for  laboratory  phantom  experiment. 

Eight  sources  and  detectors  were  used  for  light  delivery  and  detection. 

For  demonstrating  HCM  experimentally,  we  utilized  the  setup  depicted  in  Fig.  4.  A  broadband  pulsed 
laser  source  (SC-450,  Fianium  Inc.,  Eugene,  Oregon)  with  a  max  output  power  of  40  mW  (measured  at 
the  tip  of  an  optical  fiber)  was  utilized  as  the  illumination  source.  Although  the  pulsed  laser  has  a 
repetition  rate  of  20  MHz,  due  to  the  limited  frequency  response  of  the  photodiodes,  we  still  considered 
our  system  as  a  CW  system.  The  laser  output  was  collimated  and  fed  through  a  filter  wheel  where  a 
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desired  wavelength  of  780nm  was  selected.  In  order  to  calibrate  the  fluctuations  from  the  laser  source, 
a  beam  sampler  (BS)  was  utilized,  and  a  sample  beam  was  fed  into  another  photodiode  (PD)  for  further 
calibration.  The  light  was  again  coupled  into  another  optical  fiber  and  fed  into  a  multiplexor.  The 
multiplexor,  as  controlled  by  a  computer,  illuminated  8  locations  on  an  optode  probe  sequentially.  We 
employed  silica  optical  fibers  (core  diameter  of  600  urn)  for  light  delivery. 
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Figure  5:  Flow  chart  depicting  various  stages  of  detection  electronics  utilized  in  proposed 
instrumentation 

For  the  detection  electronics  (Fig.  5),  we  utilized  photodiodes  (OPT101,  Burr-Brown  Corp.,  Tucson, 
Arizona)  to  collect  the  diffused  light.  An  analog  low  pass  filter  (LPF)  was  employed  to  filter  out 
unwanted  higher  frequencies.  Next,  the  electrical  signal  was  amplified  using  an  operational  amplifier  in 
the  inverting  mode.  The  operational  amplifier’s  gain  was  controlled  by  changing  the  resistance  of 
feedback  resistor.  For  digital  control  of  gain,  a  digital  potentiometer  was  utilized  such  that  the  gain  was 
controlled  from  the  computer.  Next,  the  electrical  signal  was  further  fed  into  a  buffer  amplifier  which 
itself  was  an  operational  amplifier  in  non-inverting  mode  with  unity  gain.  Finally,  we  utilized  a  National 
Instruments  DAQ  card  with  16-bit  resolution  to  convert  analog  signals  into  digital  format  and  stored  in 
the  computer.  The  proposed  instrument  setup,  including  scanning  mechanism,  electrical  gain,  and  data 
acquisition,  was  controlled  using  Labview  software. 

Experimental  setup 

Although  in  the  computer  simulations,  we  demonstrated  our  FICM  using  a  multilayered  model,  we 
utilized  a  single  layer  phantom  for  experimental  demonstration.  We  preferred  a  single  layered  phantom 
because  of  the  complexity  involved  in  the  preparation  of  a  multilayered  phantom.  Using  a  single  layered 
phantom  with  an  absorber  inside,  we  skipped  Step  1  in  our  reconstruction  technique  and  directly 
started  the  reconstruction  from  step  2.  The  reason  for  skipping  Step  Iwas  that  it  involved  inducing  the 
prior  information;  for  this  particular  case,  we  did  not  use  any  prior  information.  The  phantom  setup  is 
depicted  in  Fig.  6.  A  homogeneous  liquid  tissue-mimicking  phantom  was  prepared  by  filling  a  container 
of  dimensions  of  15x10x10  cm3  with  1%  Intralipid  solution.  This  solution  served  as  the  homogeneous 
background  medium  with  an  absorption  coefficient  ( ya )  of  0.1  cm-1  and  reduced  scattering  coefficient 
(p’s)  of  10  cm-1.  A  spherical  absorbers  (fja  =  0.3  cm"1)  of  1-cm  diameter  were  placed  at  1.5-cm  distance 
around  the  center  of  an  optode  array  which  was  placed  from  one  side  surface  of  the  container  ,  as 
shown  in  Fig.  6. 


Figure  6:  Experimental  setup  used  in  the 
laboratory  phantom  measurement.  An  optode 
probe  has  been  placed  on  the  side  of  the  tank 
to  avoid  contact  with  the  intralipid  solution.  An 
absorber  has  been  placed  1 .5  cm  from  the  side 
of  the  tank  within  the  Intralipid  solution. 


Experimental  results 

Panels  in  Fig.  7  provide  overall  comparisons  among  the  reconstructed  images  that  were  obtained  using 
the  experimental  setup  given  in  Figs.  4-6.  Figure  7(a)  demonstrates  the  result  obtained  from  the  regular 
iterative  DOT  reconstruction  technique.  While  this  figure  shows  relatively  good  reconstruction  in  both 
location  and  size  for  the  embedded  absorber,  we  see  very  poor  recovery  in  absorption  contrast  for  the 
reconstructed  object:  the  reconstructed  value  of  iJa  was  0.012  cm'1  while  the  expected  value  was  0.03 
cm'1.  As  seen  in  Fig.  7(b),  step  2  in  our  HCM  is  able  to  determine  the  probable  location  of  the  absorber; 
Figs.  7(c)  and  7(d)  show  the  gradual  improvement  in  the  contrast  recovery.  By  the  end,  both  location 
and  absorption  contrast  of  the  embedded  absorber  were  reconstructed  much  better  to  meet  the 
expected  values. 
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Figure  7:  reconstructed  /ja  values  and  maps  in  mm'1  (a)  using  regular  iterative  DOT 
reconstruction  technique  and  HCM  after  (b)  step  2,  (c)  step  3,  and  (d)  step  4. 


3.  Key  Research  Accomplishments  and  Reportable  Outcomes 

Key  research  accomplishments: 

(1)  We  have  shown  the  feasibility  of  TRUS-DOT  through  computer  simulations  and  laboratory 
phantom  experiments. 

(2)  We  have  designed  a  TRUS-compatible,  DOT-based  image  system  (780nm),  in  which  the  photo 
diodes  were  placed  on  the  trans-rectal  probe.  Optical  signals  were  recorded  and  used  for 
estimating  the  absorption  coefficient. 

(3)  We  validated  the  system  using  laboratory  phantoms.  For  improved  image  reconstruction,  we 
have  also  developed  a  hierarchical  clustering  method  (HCM)  to  improve  the  accuracy  of  image 
reconstruction  with  limited  prior  information. 

Reportable  outcomes: 

(1)  Peter  LeBoulluec,  Hanli  Liu,  and  Baohong  Yuan,  “A  cost-efficient  frequency-domain 
photoacoustic  imaging  system,”  Am.  J.  Phys.  Vol.  81,  712  (2013);  doi:  10.1 1 19/1.4816242 

(2)  Venkaiah  C.  Kavuri  and  Hanli  Liu,  “Development  of  multispectral  transrectal  ultrasound 
compatible  near  infrared  imaging  system  for  early  detection  of  prostate  cancer,”  SPIE, 
Photonics  West,  BiOS  Biomedical  Optics  Symposium,  paper  8578-57,  Feb.  2-7,  2013,  San 
Francisco,  California. 
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4.  Conclusions 


As  mentioned  in  Section  1,  we  have  modified  our  original  goals  or  ideas  by  developing  transrectal 
ultrasound-guided  diffuse  optical  tomography  (TRUS-DOT)  for  improved  prostate  cancer  imaging.  In 
this  report,  we  have  shown  the  feasibility  of  TRUS-DOT  through  computer  simulations  and  laboratory 
phantom  experiments.  Specifically,  we  have  designed  a  TRUS-compatible,  DOT-based  image  system 
(780nm),  in  which  the  photo  diodes  were  placed  on  the  trans-rectal  probe.  Optical  signals  were 
recorded  and  used  for  estimating  the  absorption  coefficient.  We  validated  the  system  using  laboratory 
phantoms.  For  improved  image  reconstruction,  we  have  also  developed  a  hierarchical  clustering 
method  (HCM)  to  improve  the  accuracy  of  image  reconstruction  with  limited  prior  information. 

For  the  extended  period  of  funding  (till  March  31,  2014,  which  has  been  approved),  we  have  following 
modified  Statement  of  Work  (SOW)  to  complete: 

Task  1 :  We  will  make  more  tissue  phantoms  for  TRUS-compatible  imaging  experiments  by  using 
gelatin  materials  mixed  with  blood  and  intralipid.  We  will  perform  the  TRUS-DOT  imaging 
measurements  with  the  prostate-shaped  phantoms  that  simulate  cancerous  tissues. 

Task  2:  We  will  characterize  the  quality  of  reconstructed  TRUS-DOT  images,  including  image 

resolutions,  position  errors  of  the  given  simulated  tumors,  and  quantification  accuracy  of  the 
physiological  parameters  of  the  simulated  tumor  tissues. 

In  the  meantime,  we  will  publish  our  results  in  peer-reviewed  journals. 


to 
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Photoacoustic  (PA)  imaging  techniques  have  recently  attracted  much  attention  and  can  be  used  for 
noninvasive  imaging  of  biological  tissues.  Most  PA  imaging  systems  in  research  laboratories  use 
the  time  domain  method  with  expensive  nanosecond  pulsed  lasers  that  are  not  affordable  for  most 
educational  laboratories.  Using  an  intensity  modulated  light  source  to  excite  PA  signals  is  an 
alternative  technique,  known  as  the  frequency  domain  method,  with  a  much  lower  cost.  In  this 
paper,  we  describe  a  simple  frequency  domain  PA  system  and  demonstrate  its  imaging  capability. 
The  system  provides  opportunities  not  only  to  observe  PA  signals  in  tissue  phantoms  but  also  to 
acquire  hands-on  skills  in  PA  signal  detection.  It  also  provides  opportunities  to  explore  the 
underlying  mechanisms  of  the  PA  effect.  ©2013  American  Association  of  Physics  Teachers. 
[http://dx.doi.org/10.!  1 19/1.4816242] 


I.  INTRODUCTION 

Optical  imaging  of  cancer  has  been  intensively  studied  in 
recent  years  due  to  its  unique  and  high  sensitivity  to  endoge¬ 
nous  and  exogenous  tumor  contrast.  These  optical  techni¬ 
ques  are  usually  limited  either  in  penetration  depth,  e.g.,  a 
few  hundred  of  microns  for  optical  microscopy,  or  in  spatial 
resolution,  e.g.,  a  few  millimeters  for  optical  diffuse  optical 
tomography  (DOT). 13  Obviously,  a  tradeoff  exists  between 
the  imaging  depth  and  spatial  resolution.4  To  overcome  this 
limitation,  ultrasound  techniques  have  been  combined  with 
optical  approaches,  yielding  various  diagnostic  techniques 
such  as  photoacoustic  (PA)  imaging.4  Using  PA  techniques, 
the  ratio  of  imaging  depth  to  spatial  resolution  can  reach  val¬ 
ues  up  to  about  100,  which  is  ~10  times  higher  than  what 
can  be  achieved  with  conventional  DOT.  PA  techniques  can 
potentially  be  used  for  imaging  cancers  in  human  breast, 
prostate,  skin,  thyroid,  neck,  head,  and  others  areas.4  PA 
techniques  provide  optical  contrast,  which  is  usually  much 
more  sensitive  to  functional  and  molecular  information  of 
the  tissue  than  ultrasound  techniques,  making  it  possible  to 
measure  the  concentrations  of  total  hemoglobin,  oxy-  and 
deoxy-hemoglobin,  and  specific  molecules  regulating  tumor 
growth  and  metastasis.4'5  Compared  with  pure  optical  imag¬ 
ing  methods,  such  as  DOT,  PA  techniques  have  much  higher 
spatial  resolution  and  similar  imaging  depth.4'5  The  basic 


mechanism  of  the  PA  effect  can  be  briefly  explained  as  fol¬ 
lows:  (1)  tissue  absorbs  the  energy  of  light  illumination;  (2) 
the  absorbed  energy  leads  to  a  local  temperature  rise  in  the 
tissue;  and  (3)  a  pressure  or  sound  wave,  called  the  PA  sig¬ 
nal,  is  generated  due  to  the  thermoelastic  effect.4 

Most  PA  imaging  systems  in  large  research  laboratories 
use  a  nanosecond  pulsed  laser  with  high  pulse  energy.  A  short 
light  pulse  avoids  energy  loss  due  to  thennal  energy  diffusion 
before  the  generation  of  PA  signals,  which  is  known  as  ther¬ 
mal  confinement.6  A  nanosecond  pulsed  laser  with  high  pulse 
energy  typically  costs  tens  of  thousands  of  dollars  and  is  usu¬ 
ally  positioned  on  a  large  optical  table.  Low  cost  and  compact 
PA  imaging  techniques  that  do  not  rely  on  nanosecond  pulsed 
laser  systems  are  useful  for  educational  and  research  labora¬ 
tories  with  limited  resources  and  for  global  healthcare.  1 1  In 
this  study,  we  have  developed  a  frequency-domain  PA  (FD- 
PA)  imaging  system  by  using  an  inexpensive  laser  diode,  a 
conventional  single  element  ultrasound  transducer,  and  a 
lock-in  amplifier  to  improve  the  signal-to-noise  ratio. 12-16 
The  theory  and  the  experimental  implementation  of  FD-PA 
have  been  extensively  described  in  the  literature.  ’  ’  '  A 
simple,  compact,  and  cost  efficient  FD-PA  imaging  system  is 
an  excellent  teaching  and  research  tool  for  educators  and 
researchers.  It  allows  students  and  the  general  public  to 
understand  the  underlying  mechanisms  of  the  PA  effect  and 
its  applications,  such  as  cancer  imaging. 
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II.  FUNDAMENTAL  THEORY  OF  FREQUENCY- 
DOMAIN  PHOTOACOUSTICS 

The  wave  equation  used  to  quantify  photoacoustic  pres¬ 
sure  is  usually  expressed  as  follows: 


- pdH 
Cp  dt 


(1) 


where  P  represents  the  photoacoustic  signal,  H  is  the  laser 
heating  function,  v  is  the  speed  of  sound,  /i  is  the  thermal 
expansion  coefficient,  and  Cp  is  the  heat  capacity  per  unit 
mass  at  constant  pressure  of  the  irradiated  medium.  In  the 


frequency  domain,  the  laser  intensity  is  usually  modulated  as 
a  sinusoidal  function  with  frequency  /.  Assuming  the  ampli¬ 
tude  of  the  modulated  intensity  is  Iq,  the  heating  function  can 
be  expressed  as: 

H  =  nJo  exp(-icot),  (2) 

where  co  =  2nf  and  jia  is  the  absorption  coefficient  of  the 
medium,  which  is  usually  proportional  to  the  absorber  con¬ 
centration  in  the  sample.  When  considering  an  infinitely 
long  and  optically  thin  planar  or  cylindrical  absorber,  the 
generated  PA  waves  in  the  frequency-domain  can  be  repre¬ 
sented  by 
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(planar  absorber) 
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where  /  is  the  thickness  of  the  slab,  a  is  the  radius  of  the  cyl¬ 
inder,  and  r  is  the  radial  coordinate.  Jn  and  J  i  are  the  zeroth- 
and  first-order  Bessel  functions,  respectively.  //(|  and  H\ 
represent  the  zeroth-  and  first-order  Hankel  functions,  respec¬ 
tively.  q  is  the  dimensionless  frequency  and  t  is  the  dimen¬ 
sionless  time.  For  a  slab,  q  =  col/ 2v  and  t  =  2 v(t  —  ~^lJ2). 
For  a  cylinder,  q  =  coa/v  and  t  =  vt/ a.  p  =  pin/ pout  and  v 
=  Vjn  j vOUf  are  the  dimensionless  density  and  speed-of-sound 
parameters,  respectively.  The  subscripts  “in”  refers  to  inside 
the  absorber  and  “out”  refers  to  outside  the  absorber.  For  a 
small  spherical  absorber  with  radius  a,  the  PA  wave  can  be 
expressed  as 
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x  exp  (—iqt), 

where  q  =  coa/v  and  f  =  ()))  (t  —  '  v") . 


ipv  sin(^) 
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,  (cylindrical  absorber) 


(4) 


Using  Eqs.  (3) — (5),  the  following  can  be  concluded  about 
the  strength  of  the  FD-PA  signal:  (1)  it  is  a  sinusoidal  func¬ 
tion  of  time  with  the  same  frequency  co  as  that  of  the  modu¬ 
lated  light;  (2)  it  is  proportional  to  the  optical  absorption 
coefficient  pa  of  the  medium;  (3)  it  is  proportional  to  the  am¬ 
plitude  of  the  modulated  intensity  Iq\  and  (4)  it  depends  on 
the  modulation  frequency  co  via  <7.IN  ly  In  this  study,  the  sam¬ 
ple  setup  can  be  assumed  to  be  a  cylindrical  absorber, 
described  by  Eq.  (4). 

III.  MEASUREMENT  SYSTEM  AND  DATA 
PROCESSING 

A.  Measurement  system 

Figure  1  shows  the  experimental  setup.  A  function  genera¬ 
tor  (FG,  Agilent  33120A,  Agilent  Tech)  generates  a  sinusoi¬ 
dal  voltage  signal  of  frequency  /  and  a  synchronized  TTL 
(Transistor-Transistor  Logic)  signal  with  the  same  frequency 


Fig.  1 .  The  experimental  setup  used  for  PA  imaging.  The  tank  with  the  ICG  sample  is  filled  with  either  water  or  an  intralipid  solution.  The  function  generator 
provides  a  sinusoidal  voltage  signal  to  the  driver  of  the  laser  diode.  A  synchronized  TTL  signal  is  sent  to  the  lock-in  amplifier  and  serves  as  the  reference  sig¬ 
nal.  The  laser  diode  emits  a  modulated  light  beam  that  is  focused  on  the  sample.  A  shutter  is  used  to  block  the  laser  in  order  to  measure  the  background  noise. 
The  modulated  light  is  absorbed  by  the  sample,  causing  the  emission  of  PA  waves,  which  are  detected  by  the  ultrasound  transducer.  The  output  voltage  signal 
from  the  transducer  is  amplified,  filtered,  and  sent  to  the  lock-in  amplifier. 
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Fig.  2.  Schematic  of  the  modulation  circuit  used  in  our  setup.  The  function  generator  provides  DC  and  AC  power  to  the  circuit.  The  resistor  and  the  potentiometer 
are  chosen  to  provide  a  modulated  current  within  the  current  range  of  the  laser  diode.  The  photodiode  detector  and  oscilloscope  are  used  to  verify  the  modulation. 


and  a  fixed  phase  shift.  The  sinusoidal  signal  is  sent  to  a 
homemade  circuit,  shown  in  Fig.  2,  to  drive  the  laser  diode 
(L785P100,  Thorlabs).  Thus,  the  laser  intensity  is  modulated 
with  the  frequency  /.  The  laser  diode  has  a  central  wave¬ 
length  of  785  nm  and  a  power  of  ~100  mW  when  operating 
in  DC  mode.  A  lens  is  used  to  collimate  the  laser  beam  and  a 
mechanical  shutter  is  used  to  manually  block  the  laser  illu¬ 
mination  to  measure  the  background  signal.  The  synchron¬ 
ized  TTL  signal  from  the  FG  is  sent  to  a  lock-in  amplifier 
(LIA,  SR844,  Stanford  Research  Systems)  and  used  as  a  ref¬ 
erence  signal.  The  sample  and  ultrasound  transducer  (UST, 
Olympus  NDT,  1  in.  inch  focal  length)  are  submerged  into  a 
transparent  10-gal  tank  (20  cm  wide,  40  cm  long,  and  24  cm 
high),  filled  with  either  water  or  an  intralipid  solution.  The 
sample  used  is  an  indocyanine  green  (ICG)  aqueous  solution, 
which  is  injected  into  a  partially  optically  and  acoustically 
transparent  tube  (MRE-095,  Braintree  Scientific)  that  is  ver¬ 
tically  positioned  in  the  tank.  The  outer  and  inner  diameters 
of  the  tube  are  2.5  and  1.7  mm,  respectively.  The  UST  is 
focused  on  the  sample.  Both  the  lateral  and  axial  sizes  of  the 
focal  zone,  normally  measured  in  terms  of  the  full  width  at 
half  maximum  (FWHM),  depend  on  the  central  frequency  of 
the  UST.  When  the  intensity-modulated  laser  reaches  the 
optically  absorbing  sample,  a  PA  wave  is  generated.  The  PA 
wave  is  detected  and  converted  into  a  voltage  signal  by  the 
UST  and  further  amplified  (ZFL-1000LN,  Mini  Circuits)  and 
filtered  (SLP-5+,  Mini  Circuits).  A  total  gain  of  40  dB  is 
applied  by  using  two  identical  amplifiers  in  series,  driven  by 
a  DC  power  supply  (BK  Precision  1506).  The  low  pass  filter 
(SLP-5+,  Mini  Circuits)  is  not  required  because  the  LIA  has 
a  very  narrow  bandwidth.  The  processed  PA  signal  is  deliv¬ 
ered  to  the  LIA  and  the  amplitude  of  the  PA  signal  and  the 
phase  difference  between  the  PA  signal  and  the  reference 
signal  are  displayed  on  the  screen  of  the  LIA. 

Figure  2  shows  the  principle  of  the  modulation  circuit.  By 
generating  an  AC  signal  with  a  DC  offset,  the  FG  can  serve 
as  a  DC  and  AC  power  source.  A  resistor  (Rl)  and  a  potenti¬ 
ometer  are  used  to  limit  and  control  the  current  flowing  into 
the  laser  diode,  respectively.  The  laser  diode  has  a  typical 
threshold  current  of  35  mA.  The  DC  offset  of  the  FG  is 
Vqc  =  2  V  and  the  AC  peak-to-peak  is  Vac  =  4  V.  The  total 
resistance  of  Rl  and  the  potentiometer  is  about  1 1  Q.  A  pho¬ 
todiode  (EOT,  ET-2030A)  is  used  to  verify  the  modulation 
of  the  laser  intensity  via  an  oscilloscope  (2530B  Digital 
Storage  Oscilloscope,  BK  Precision). 

B.  Data  processing 

It  is  common  that  the  LIA  shows  a  background  signal  even 
when  the  laser  is  turned  off  or  blocked.  This  background  signal 
is  mainly  caused  by  electronic  interference  from  the  driving 
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signal  generated  by  the  FG.  Fortunately,  this  background  sig¬ 
nal  is  independent  of  the  location  of  the  UST  or  the  sample.  A 
simple  way  to  eliminate  the  effect  of  this  background  signal  is 
to  mathematically  subtract  it  from  the  measured  PA  signal. 
Figure  3  schematically  shows  the  relationship  between  the 
three  signals,  based  on  the  amplitude  and  the  phase  measured 
with  the  LIA:  (1)  the  measured  signal,  (2)  the  background  sig¬ 
nal,  and  (3)  the  PA  signal.  The  angle  between  the  vector  repre¬ 
senting  each  signal  and  the  horizontal  axis  is  the  phase.  The 
PA  signal,  vector  (3)  in  Fig.  3,  is  the  difference  between  the 
measured  signal  and  the  background  signal.  The  magnitude  of 
the  PA  signal  is  correlated  with  the  optical  absorption  coeffi¬ 
cient  of  the  sample.  All  the  calculations,  image  processing,  and 
plotting  were  carried  out  using  matlab  (Mathworks),  which  is 
commonly  used  in  engineering  laboratories.  Other  software 
packages,  such  as  mathcad,  can  be  used  as  alternatives. 

IV.  RESULTS  AND  DISCUSSIONS 

Before  the  PA  measurements,  the  laser  beam  was  posi¬ 
tioned  to  be  coaxial  with  the  UST.  The  tank  was  filled  with  a 
1%  intralipid  solution  to  simulate  biological  tissues.  The 
absorption  and  scattering  coefficients  of  the  intralipid  solution 
are  0.04  cm”1  and  8.4  cm-1,  respectively,  as  measured  with  an 
ISS  Oximeter.  The  tube  was  filled  with  an  ICG  aqueous  solu¬ 
tion,  with  a  concentration  of  0.5  g/1,  to  simulate  an  absorbing 
target.  ICG  is  a  relatively  weak  fluorophore  in  the  near  infra¬ 
red  range  with  a  quantum  yield  of  <~1%  in  aqueous  media20 
and  has  been  used  as  an  optical  absorber  for  PA  imaging."1 
The  laser  diode  and  the  UST  were  mounted  on  the  same  trans¬ 
lation  stage  and  their  relative  position  was  thus  fixed.  By  mov¬ 
ing  the  translation  stage,  the  PA  signal  distributions  along  the 
x  axis  and  in  the  x-y  plane  were  acquired.  Although  the  co¬ 
axial  setup  between  the  laser  beam  and  the  UST  is  not 
required,  it  improves  the  signal-to-noise  ratio. 


Fig.  3.  The  relation  between  the  measured  signal  (vector  1),  the  background 
(vector  2),  and  PA  (vector  3)  signals. 
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Fig.  4.  The  normalized  PA  signal  strength  as  a  function  of  the  horizontal  position  of  the  sample  relative  to  the  co-axis  of  the  laser  beam  and  the  UST.  The  data 
were  obtained  with  a  1  MHz  (dashed  curve)  and  a  2.25  MHz  (solid  curve)  UST.  The  dashed  vertical  lines  show  the  inner  diameters  and  the  solid  vertical  lines 
show  the  outer  diameters  of  the  phantom  tube  that  is  filled  with  0.5  g/1  ICG  solution  at  a  depth  of  1.9  mm  in  a  1%  intralipid  solution.  The  diameter  of  the  colli¬ 
mated  laser  beam  is  ~  1.0  mm.  The  error  bars  represent  the  standard  deviation  of  4  measurements. 


Figure  4  shows  the  strength  of  measured  PA  signal  as  a 
function  of  the  horizontal  position  of  the  tube  relative  to  the 
axis  of  the  laser  beam  and  the  UST.  The  dashed  and  solid 
vertical  lines  represent  the  inner  and  outer  diameter  of  the 
tube,  respectively.  Clearly,  the  strength  of  the  PA  signals  rise 
when  the  laser  beam  gradually  moves  into  the  tube  region 
and  fall  when  the  laser  beam  moves  away  from  the  tube 
region.  These  results  indicate  that  the  highly  absorbing  ICG 
tube  generates  significant  PA  signals  compared  to  the  sur¬ 
rounding  intralipid  solution  that  has  a  much  lower  absorption 
coefficient.  The  FWHM  of  the  PA  data  using  the  1  MHz 
UST,  shown  in  Fig.  4,  is  about  2.2  mm.  It  is  larger  than  the 
inner  diameter  of  the  tube  (1.7  mm)  and  the  UST  lateral  size 
( 1 .3  mm)  since  the  FWHM  is  mainly  determined  by  the  con¬ 
volution  of  the  profiles  of  the  cross  section  of  the  tube  and 
the  lateral  focal  zone  of  the  UST. 

When  the  frequency  of  the  UST  is  increased  from  1  to 
2.25  MHz,  the  FWHM  of  the  UST’s  lateral  focal  zone  is 
reduced  from  ~1.3  to  ~0.72mm.  The  PA  spatial  resolution 
can  be  thus  improved  by  using  a  higher  frequency  UST. 
Figure  4  shows  that  the  FWHM  of  the  PA  signal  measured 
with  the  2.25  MHz  UST  is  0.98  mm  which  is  smaller  than  the 
2.2  mm  width  measured  with  the  1  MHz  UST.  However,  the 


Fig.  5.  The  maximum  PA  signal  strength  as  a  function  of  the  tube  depth 
acquired  with  a  2.25  MHz  UST  in  a  1%  intralipid  solution.  The  diameter  of 
the  laser  beam  is  ~  1.0  mm  and  the  ICG  concentration  is  0.5  g/1-  The  error 
bars  represent  the  standard  deviation  of  4  measurements. 


width  is  smaller  than  the  tube  size  and  the  reason  is  unclear. 
One  possible  reason  may  be  that  the  PA  signal  generated 
from  regions  close  to  the  edge  of  the  tube  is  too  weak  to  be 
detected.  Higher  frequencies  provide  smaller  focal  sizes  and 
higher  spatial  resolution.  A  UST  with  a  central  frequency 
between  1  and  10  MHz  achieves  an  acceptable  spatial  resolu¬ 
tion.  At  frequencies  below  1  MHz,  the  resolution  will  be 
degraded.  Frequencies  above  10  MHz  will  also  work  but 
increases  the  cost  and  the  complexity  of  the  system. 

When  the  distance  between  the  left-side  of  the  tank  wall 
and  the  left  outer  surface  of  the  tube,  defined  as  the  depth  of 
the  tube,  increases,  the  strength  of  the  PA  signal  decreases. 
This  is  a  consequence  of  the  reduction  in  the  light  intensity 
reaching  the  sample  due  to  the  scattering  in  the  medium.  The 
dependence  of  the  strength  of  the  PA  signal  on  the  thickness 
or  depth  of  the  medium  is  shown  in  Fig.  5. 

Increasing  the  ICG  concentration  increases  the  absorption 
coefficient  of  the  tube,  which  raises  the  PA  signal  strength. 
Figure  6  shows  the  peak  strength  of  the  PA  signal  as  a  func¬ 
tion  of  the  ICG  concentration.  The  PA  signal  strength 
increases  when  the  concentration  is  increased  but  appears  to 
saturate  at  high  ICG  concentrations.  This  may  indicate  that 
the  light  cannot  penetrate  into  a  highly  concentrated  ICG 


Fig.  6.  The  maximum  PA  signal  strength  as  a  function  of  the  ICG  concentra¬ 
tion  acquired  with  a  2.25  MHz  UST  in  a  1%  intralipid  solution  at  a  depth  of 
1.9  mm.  The  diameter  of  the  collimated  laser  beam  is  ~  1.0  mm.  The  error 
bars  represent  the  standard  deviation  of  5  measurements. 
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Fig.  7.  (a)  A  C-mode  ultrasound  image  and  (b)  a  C-mode  PA  image  of  the  ICG  filled  tube.  The  UST  frequency  is  2.25  MHz,  the  ICG  concentration  is  2g/l, 
and  the  intralipid  concentration  is  1%.  (c).  An  ultrasound  B-mode  image  showing  the  cross  section  of  the  same  ICG  filled  tube.  The  white  dotted  line  indicates 
the  depth  of  the  C-mode  images  in  (a)  and  (b).  Note  that  the  depth  of  zero  in  this  figure  is  the  surface  of  the  UST,  which  is  different  from  the  definition  of  the 
tube  depth),  (d)  A  diagram  to  show  the  B-mode  and  C-mode  imaging  planes.  The  B-mode  image  is  parallel  to  the  x-z  plane  and  the  C-mode  is  parallel  to  the 
x-y  plane  (close  to  the  right-side  edge  of  the  tube). 


solution  due  to  the  large  absorption  coefficient,  suggesting 
that  the  light  absorption  is  confined  to  a  limited  region  of  the 
ICG  solution. 

Figures  7(a)  and  7(b)  show  two-dimensional  (2D)  ultra¬ 
sound  and  PA  images  of  the  tube  in  the  x-y  plane  (C-mode, 
see  Fig.  7(d)  for  the  coordinates  and  the  imaging  plane).  The 
ultrasound  image  was  obtained  with  the  same  UST  using  the 
conventional  pulse-echo  method.  The  UST  was  connected  to 
a  pulse-generator-receiver  (5073PR  pulser/receiver,  Olympus) 
and  the  reflected  ultrasound  data  was  acquired  with  the  oscil¬ 
loscope.  The  UST  was  raster  scanned  in  the  x-y  plane.  The 
laser  was  turned  off  because  it  is  not  used  in  ultrasound  imag¬ 
ing.  The  reflected  ultrasound  echo  signal  at  each  location,  also 
called  the  A-line,  included  the  acoustic  information  along  the 
z  direction.  After  scanning  the  x-y  plane,  all  A-lines  can  be 
used  to  construct  a  3D  image  of  the  tube.  Figures  7(a)  shows  a 
2D  image  in  the  x-y  plane,  the  so-called  a  C-mode  ultrasound 
image,  at  a  depth  indicated  by  the  white  line  in  Fig.  7(c). 
Figure  7(c)  shows  a  2D  image  in  the  x-z  plane,  the  so-called  a 
B-mode  ultrasound  image,  representing  the  cross  section  of 
the  tube.  The  PA  image  shown  in  Fig.  7(b)  was  obtained  with 
the  system  described  in  Fig.  1 .  This  system  does  not  have  the 
capability  to  resolve  the  PA  signal  along  z  direction  because 
the  laser  diode  is  continuously  modulated.  This  is  one  disad¬ 
vantage  of  the  PA  technique  compared  with  the  time-domain 
method.  Therefore,  only  a  2D  image  in  the  x-y  plane  is  shown 
in  Fig.  7(b). 
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The  data  points  shown  in  Figs.  7(a)  and  7(b)  were  obtained 
by  raster  scanning  the  UST/laser  system  in  the  x-y  plane. 
The  step  sizes  in  the  lateral  (x)  direction  and  in  the  vertical 
(z)  direction  were  0.5  mm  and  1.25  mm,  respectively.  Each 
rectangle  represents  the  intensity  of  a  single  measurement. 
The  two  solid  and  dashed  white  vertical  lines  represent  the 
outer  and  inner  diameter  of  the  tube,  respectively.  The  aver¬ 
age  FWHM  of  the  ultrasound  image  shown  in  Fig.  7(a)  is 
about  0.98  mm,  which  is  smaller  than  the  inner  diameter  of 
the  tube.  Figure  7(b)  shows  the  corresponding  PA  image  in 
the  x-y  plane.  The  image  clearly  shows  the  tube  with  optical 
(absorption)  contrast.  The  average  FWHM  is  about  0.88  mm, 
which  is  also  smaller  than  the  inner  diameter  of  the  tube  for 
reasons  that  are  not  clear. 

Figure  7(c)  shows  the  ultrasound  image  of  the  cross- 
section  of  the  tube.  To  be  consistent  with  conventional  ultra¬ 
sound  imaging  rules,  zero  depth  in  Fig.  7(c)  is  defined  as  the 
surface  of  the  UST  because  that  is  the  location  where  the 
ultrasound  pulses  are  generated;  note  that  this  is  different 
from  the  definition  of  the  tube  depth.  The  two  bright  areas 
show  the  two  boundaries  of  the  tube.  Note  that  the  inner 
boundaries  of  the  tube  cannot  be  resolved  from  this  image 
due  to  the  limited  resolution  of  the  UST  (2.25  MHz).  The 
dotted  horizontal  line  in  Fig.  7(c)  indicates  the  depth  of  the 
ultrasound  C-mode  image  shown  in  Fig.  7(a).  Clearly,  the 
ultrasound  image  shows  the  boundaries  between  the  tube  and 
the  surrounding  medium. 
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Table  I.  The  major  components  of  the  PA  system  and  their  features,  costs,  and  possible  substitutes.  The  listed  prices  are  for  guidance  only. 


Item 

Model 

Feature 

Cost  (USD) 

Substitutes 

Lock-in  amplifier 

SRS  844 

Broadband  200  MHz 

~$8,000 

Regular  LIA,  <$2000 

Function  generator 

Agilent  33120  A 

Multiple  function  and 

20  MHz  bandwidth 

~$1600 

Regular  FG,  <$1000 

Focused  ultrasound 

transducer 

Olympus-NDT  V314-SU-F 

Point  focused,  NA=0.375 

~$350 

Other  focused  USTs  with  a  central  frequency 
between  1  and  10  MHz  and  reasonable  sensitivity 

Laser  diode 

Thorlabs  L785P100 

785  nm  and  100  mw  CW 

~$40 

Other  laser  diodes  with  power  ~  100  mw 

RF  amplifiers  (2x) 

Minicircuits  ZFL-1000LN 

Wideband,  0.1 -1000  MHz 

—$178 

(2x  at  $89  each) 

Other  RF  amplifiers  with  a  gain  ~40  dB 

Low  pass  filter 
(not  required) 

Minicircuits  SLP-5+ 

Cutoff  at  5  MHz 

~$35 

Other  low  pass  filters  with  an 
appropriate  cutoff  frequency 

V.  COST  AND  COMPACTNESS  ANALYSES 

Table  I  shows  the  general  costs,  features,  and  possible 
substitutes  of  the  major  components  of  the  current  PA  sys¬ 
tem.  The  cost  of  the  system  is  dominated  by  the  cost  of  the 
LIA  and  the  FG  although  lower  cost  LIAs  and  FGs  could  be 
adopted.  Clearly,  the  cost  of  a  FD-PA  imaging  system  is 
much  less  expensive  than  that  of  a  time-domain  imaging  sys¬ 
tem.  The  cost  of  the  matlab  software  is  not  included  in 
Table  I.  A  student  version  of  matlab  is  relatively  cheap  but 
other  free  or  relatively  inexpensive  software  packages  can  be 
used  as  alternatives.  All  components  of  the  FD-PA  system 
are  compact.  Although  we  have  not  attempted  to  reduce  the 
overall  size  of  our  setup  in  the  current  study,  it  is  possible  to 
integrate  the  entire  system  into  a  small  and  portable  box  by 
replacing  the  FG  and  LIA  with  a  small,  customized  circuit 
board.  It  will  be  extremely  difficult  for  a  time-domain  PA 
system  to  achieve  the  same  compactness. 

VI.  CONCLUSIONS 

A  FD-PA  imaging  system  was  designed  and  its  imaging 
capability  was  demonstrated  with  tissue-like  phantoms. 
The  imaging  principle  and  data  processing  method  were 
discussed.  The  measured  PA  signal  strength  depends  on 
the  absorption  coefficient,  the  depth  of  the  target,  and  the 
modulated  light  intensity  of  the  laser  diode.  Appropriate 
processing  of  the  measured  PA  data  and  background  inter¬ 
ference  is  an  important  key  to  correctly  display  the  PA 
images.  The  lateral  spatial  resolution  of  PA  images  is  de¬ 
pendent  on  the  ultrasound  frequencies  and  focal  sizes  of 
the  ultrasound  transducers.  Although  the  current  system 
has  a  low  axial  resolution  due  to  the  adopted  ultrasound 
transducer,  it  can  be  significantly  improved  by  using  a 
high  numerical  aperture  ultrasound  transducer  or  by  adopt¬ 
ing  a  frequency-swept  technique.9  The  FD-PA  imaging 
system  is  cost  effective  compared  to  a  time-domain  imag¬ 
ing  system. 
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